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ABSTRACT 



Context. 

Aims. We present a deep optical/near-infrared imaging survey of the Serpens molecular cloud. This survey constitutes the comple- 
mentary optical data to the Spitzer "Core To Disk" (c2d) Legacy survey in this cloud. 

Methods. The survey was conducted using the Wide Field Camera at the Isaac Newton Telescope. About 0.96 square degrees were 
imaged in the R and Z filters, covering the entire region where most of the young stellar objects identified by the c2d survey are 
located. 26524 point-like sources were detected in both R and Z bands down to R »24.5 mag and Z «23 mag with a signal-to-noise 
ratio better than 3. The 95% completeness limit of our catalog corresponds to 0.04 M Q for members of the Serpens star forming region 
(age 2 Myr and distance 260 pc) in the absence of extinction. Adopting the typical extinction of the observed area (Ay ~7 mag), we 
estimate a 95% completeness level down to M»0. 1 M Q . The astrometric accuracy of our catalog is 0.4 arcsec with respect to the 
2MASS catalog. 

Results. Our final catalog contains J2000 celestial coordinates, magnitudes in the R and Z bands calibrated to the SDSS photometric 
system and, where possible, JHK$ magnitudes from 2MASS for sources in 0.96 square degrees in the direction of Serpens. This 
data product has been already used within the frame of the c2d Spitzer Legacy Project analysis in Serpens to study the star/disk 
formation and evolution in this cloud; here we use it to obtain new indications of the disk-less population in Serpens. 
Conclusions. 

Key words, stars: catalogs - stars: formation - stars: low-mass, brown dwarfs - ISM: clouds - ISM: individual objects: Serpens 



1. Introduction 

The Spitzer Legacy Survey "Mo lecular Cores to Planet Forming 
Disks" (c2d; lEvans et all I2003I) offered a singular opportunity 
for a major advance in the study of star and planet formation. 
Thanks to its sensitivity and vawelength coverage, Spitzer al- 
lowed us to address for the first time long-standing challenges 
such as disk formation and dispersal, the physical and chemi- 
cal evolution of the circumstellar material and, in particular, to 
probe the inner planet-forming region o f disks on the ba sis of 
statistically significant samples (see, e.g jLada et al.ll2006h . 

One of the star-forming regions included in the c2d survey is 
the Serpens molecu lar cloud. Because of its prox imity (260 pc; 
Strai zvs et alJ [1996) and young age (2-6 Myr; Oliveir a et al.1 



2009>; this cloud is particularly well-suited for studies of very 



young low-mass stars and sub-stellar objects. The c2d survey in 
Serpens has provided evidence of sequential star formation in 
this cloud progressing from SW to NE and culmina ting in the 
main Serpens Core with i ts cluster of Class objects (Kaa s et al.l 
l2004t lHarvev et"ai1 120071) : moreover, the surface density of 
young stars in this region is much higher, by a factor of 10-100, 
than t hat of the other st ar-forming regions mapped by the c2d 
team (Ev ans et~aT1 12009) and includes 22% of the c2d sources 
classified as "transitional" disks. This makes Serpens the best 
region for obtaining a complete, well-defined sample of multi- 



Send offprint requests to: L. Spezzi, e-mail: lspezzi@rssd . esa . int 



wavelength observations of young stars and sub-stellar objects 
in a possible evolutionary sequence to build up a "template" 
sample for the study of disk evolution up to a few Myr within 
a single, small, well defined region. To this aim, the c2d Team 
has conducted several surveys, from X-ray to millimeter wave- 
lengths, and spectroscopic follow-ups of the newly discovered 
population of young stars in Serpens, making this cloud only the 
third star-forming region af ter Taurus and IC 348 for which such 
an unbiased dataset exists ( Goodmanl|2004t lHarvev et aLll2007t 
lEnoch et alj|2007tl01iveira et al.ll2009M20iri . 

In this paper we present the optical/near-infrared (NIR) 
imaging data collected within the frame of the c2d survey in 
Serpens. These data represent one of the critical ingredients to 
anchor the studies of envelopes and disks to the properties of 
the central stars, including those without infrared excess. Indeed, 
deep XMM-Newton data of the same field revealed a sample of 
new sources, mostly candidate weak-line T Tauri stars (WTTSs), 
half o f which hav e no counterpart in oth er catalogs dBrown et al.l 
2010). Moreover, Comeronet al. (2009) recently reported on a 
large-scale optical survey of the Lupus star-forming complex 
which, in combination with NIR data from the 2MASS cata- 
log, unveiled a large population of stars and brown dwarfs (BDs) 
that have lost their inner disks on a timescale of a few Myrs or 
less. This discovery stresses the important unknowns that still 
persist in the observational characterisation of young very low- 
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mass objects and in the time-scales and mechanisms for disk 
dissipation. 

The outline of this paper is as follows: in Sect.[2]and[3]we de- 
scribe the observations and data reduction procedure, with par- 
ticular emphasis on the photometric completeness and astromet- 
ric accuracy of the final optical/NIR catalog. In Sect. |4]we dis- 
cuss the use of this catalog within the frame of the c2d Team's 
analysis in Serpens and investigate its disk-less population. Our 
conclusions are drawn in Sect. [5] 

2. Observations 

The imaging observations presented in this work were carried 
out from 13 to 17 May 2008 using the Wide Field Camera 
(WFC) at the 2.5m Isaac Newton Telescope (INT), which is lo- 
cated at the Roque de Los Muchachos Observatory (La Palma, 
Spain). The data were collected by Ignas Snellen as part of an 
observing practicum for Leiden University astronomy students. 

The WFC is a four-chip mosaic of thinned AR coated EEV 
4Kx2K devices. Each CCD has an useful imaging area of 
2048x4 100 pixels with a pixel scale of 0.333 arcsec/pixel, cover- 
ing a total field of view 34x34 arcmin with small gaps of ~ 20" 
between adjacent chips. The average CCD read-out noise and 
gain are 7 e" and 2.8 e~/ADU, respectively. 

A total of 3 fields were observed in the direction of Serpens 
(R.A.=18' ! 29 m 49 I , Dec.=+01' / 14" 1 48 s ), covering the entire re- 
gion where most of the young stellar objects (YSOs) identi- 
fied by the Spitzer c2d survey are located (Fig. [TJ. Images 
were obtained through the R (/l c =6260 A, FWHM=1380 A) 
and Z (,l e =9100 A, FWHM=1370 A) filters approximately cov- 
ering the Sloan Dig ital Sky Survey (SDSS) r' and z' bands 
(Fukugita et al.lll996h . The total sky-area observed in each filter 
is 0.96 square degrees. Each observation in each filter was split 
into several individual exposures (ditherings), shifting the tele- 
scope pointing by ~1 arcmin between consecutive exposures; 
this allows us to cover the gaps between the CCDs and avoid 
saturation of bright sources in the field. The exposure time for 
each dithering was 120 sec in both the R and Z bands. In or- 
der to recover the photometry for bright sources saturated in 
these long-time exposures, a series of short-time exposures (i.e. 
10 sec per dithering in both filters) of the same sky-area was 
also performed. The summary of these observations is reported 
in Table Q] A standard Landolt field was also observed in both 
filters for absolute flux calibration purposes (Sect. l3~3l . 

3. Data reduction 

The raw images were processed using the IRAlfl mscred pack- 
age and a number of scripts ad hoc developed both under IRAF 
and under ID lfl We followed the s teps for the WFC data reduc- 
tion pipeline (llrwin & Lewisl2 001) developed by the Cambridge 
Astronomical Survey Unit (CASU), responsable for the process- 
ing an d archiving of the datas et from the INT Wide Field Survey 
(WFS: lMcMahon et al.ll200lh . 

3.1. Pre-reduction 

Since the bias level tends to vary, the images were first corrected 
for overscan and trimmed using the bias and trim sections as 

1 IRAF is distributed by NOAO, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under contract to the 
National Science Foundation. 

2 Interactive Data Language. 



Table 1. Journal of the observations. 



Field 


Date 


Filter 


Tcxp 


Seeing 


Air Mass 


(R.A.,Dec.) 


(d/m/y) 




(min) 


(") 




Serpens. 1 


13/05/08 


R 


120x4 


1.5 


1.50 


(18:28:43,+00:24:46) 


17/05/08 


R 


10x1 


2.0 


1.14 




15/05/08 


Z 


120x9 


1.5 


1.17 




17/05/08 


Z 


10x1 


2.1 


1.14 


Serpens_2 


14/05/08 


R 


120x9 


1.6 


1.15 


(18:28:57,+01:02:54) 


17/05/08 


R 


10x1 


2.0 


1.15 




15/05/08 


Z 


120x4 


1.5 


1.16 




17/05/08 


Z 


10x1 


1.5 


1.15 


Serpens_3 


14/05/08 


R 


120x9 


1.6 


1.15 


(18:28:36,-00:09:51) 


17/05/08 


R 


10x1 


1.6 


1.15 




17/05/08 


Z 


10x1 


2.3 


1.15 


L 104+ 


13/05/08 


R 


2x1 


1.2 


1.22 


(12:41:68,-00:34:11) 


13/05/08 


Z 


5x1 


1.3 


1.19 



1 Standard star field. 



specified in the FITS headers. For each night in which observa- 
tions for our program were performed, twilight flat frames were 
then combined to obtain the night master flat, which was then 
used to correct the science images. The corrected images are es- 
sentially linear to ~2% over the full range. Bad pixels and partial 
columns have been replaced using the bad pixel mask files avail- 
able from the WFS homepag^l 

The Z-band images suffer from significant sky fringes. In 
order to remove them, we subtracted from each Z-band image 
the fringing pattern frame available from the WFS homepage 
scaled by a specific factor to account for the amplitude of the 
fringes in the individual science frames. 

3.2. Astrometry and co-addition of images 

The astrometric calibration and relative flux scaling be- 
tween ditherings were obtained using the c-version of 
ASTROMETRDfl (M. Radovich, private communication). This 
tool performs a global astrometric solution that takes overlap- 
ping sources falling on adjacent CCDs in different ditherings 
into account. For each pointing, the astro metric solution was 
computed using the USNO-B1.0 catalogue dMonet et al.ll2003l) 
as a reference. Within the global astrometry process, the astro- 
metric solution was constrained for each CCD by both the po- 
sitions from the USNO-B1.0 catalogue and those from overlap- 
ping sources in all the other CCDs. 

The co-addition of the dithered images for a g iven filter an d 
pointing was performed using the SWARP tool (lBertinll2008l) . 
The final stacked image is a 6k x 6k frame where each pixel 
value is the median flux of the co-added ditherings normalised 
to the total exposure time and relative to the airmass and atmo- 
spheric transparency of the first frame in the dithering set. The 
absolute astrometric precision of our images is about 0.4 arcsec, 
slightly lower than RMS accuracy of the USNO-B1.0 catalogue 
(0.2 arcsec); the astrometric precision has been also confirmed 
by a cross-check with the 2MASS point-source catalog (Fig. [2j. 
The internal RMS, computed from overlapping sources in differ- 
ent exposures, is within 0.05 arcsec, indicating the good perfor- 
mance of ASTROMETRIX. 



3 http://www.ast.cam.ac.uk/ wfcsur/technical/pipeline/ 

4 See also: http://www.na.astro.it/~radovich/wifix.htm 
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Fig. 1. IRAS 100/im dust emission map of the Serpens dark cloud. The contours, from 55 to 190 MJy x sr~' in steps of 9 MJy x sr~' , 
are also drawn (white lines). The continuous line defines the area covered by the three INT+WFC pointings, while the dashed line 
is the area observed by Spitzer as part of the c2d Survey. The stars are the YSO candidates identified by the Spitzer c2d Survey. 



3.3. Photometric calibration 

Instrumental magnitu des were reported to the standard SDSS 
photometric system (Fukug ita et alJ 1 19961). To this aim, the 
Stetson's standard star field L 104 dStetson] 120001) was ob- 
served in the RZ filters. By using the IRAF package photcal, we 
first performed the aperture photometry for the standard stars, 
obtaining their instrumental magnitudes (ro and zo) corrected 
for atmospheric extinction and normalised to the exposure time. 
Then, the transformation coefficients, namely zero point (ZP) 
and colour term (c), from the WFC-INT system to the SDSS 
standard system were determined by a linear fitting of the fol- 
lowing equations: 

r' = r + c R ■ (ro -zo) + ZP R (1) 

z' =zo + c z -(ro-zo)+ZPz (2) 

where r' and z' are the standard magnitudes of Landolt's stars in 
the SDSS system. 

The mean transformation coefficients determined in our ob- 
serving run are reported in Table [2] 

Using these coefficients, instrumental magnitudes (see 
Sect, \3.4i for the observed sources in Serpens have been con- 
verted to the SDSS photometric system, which is nearly an AB 



system. Thus, magnitudes in our catalog can be turned into flux 
densities using the correction from SDSS zeropoints to AB ze- 
ropoints and the AB zeropoint flux densitjU 

correction ZPr = mag 

F° r , = 2.7769E-12 W ■ ctn 2 ■ //m -1 

correction ZPz - 0.02 mag 

F° = 1.3153E-12 W ■ cm- 2 ■ pin- 1 

3.4. The catalog extraction 

The source extraction and photometry from each stacked image 
in each filter wer e performed by using th e 2.5 version of the 
SExtractor tool bv lBertin & Arnoutsl (Il996l) . SExtractor exploits 
the aperture photometry technique, which is the faster and best 
approach for uncrowded fields such as the one considered here. 
In particular, we adopted the SExtractor adaptive aperture mag- 
nitudes (magouto), which are estimated from a flexible elliptical 
aperture around each detected object; this method is expected to 
give the most precise magnitudes for very faint objects which 

5 See http://www. sdss. org/dr6/algorithms/fluxcal. html. 
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Fig. 2. Residuals of the coordinates obtained by us for the 
sources in the surveyed area in Serpens with respect to those 
from the 2MASS catalog. 

Table 2. Mean photometric calibration coefficients for our 
WFC-INT observing run. 



Filter K* ZP 



R 0.0734 25.30±0.03 -0.033±0.024 
Z 0.0103 23.52±0.03 0.021±0.022 



i.or 




o.o - 



14 16 18 20 22 24 26 

R 

1.0 p — ' — ' — ' — 1 — ' — ' — 1 — 1 — ' — ' — ' — 1 — ' — ' — ' — 1 — ' . 1 

o.8^ 

- Amag = exp[0.934 X (mog-23.895)] + 0.0286989 ■ . f ' ; ' - 
0.6- '■■/ - 




14 16 18 20 22 24 

Z 

Fig. 3. Photometric errors versus magnitudes and relative ex- 
ponential fits for all the point-like sources detected in the 
0.96 square degrees area surveyed in Serpens. Photometry is 
from our long-time exposures (Sectj2j. 

was further complemented with near-infrared JHK$ photome- 
try by cross-matching it with t he 2MASS Point Sour ce Catalog 
(limiting magnitude K$ ~15.5; Skrutsk ie et al.l l2006). A match- 
ing radius of 2 arcsec was defined, which appears to be sufficient 
given that more than 90% of the matches corresponded to differ- 
ences between the 2MASS positions and those determined by us 
of less than 1 arcsec. When more than one 2MASS source was 
found within the 2 arcsec circle, the one closer to the position 
derived from our observations was chosen as the counterpart. 



f K = mean atmospheric extinction coefficients for La Palma. 

* ZP = photometric zero point (see Eq.[2). 

* c = photometric colour term (see Eq.|2). 



may appear extended. The detection threshold was set in order 
to select all the sources having a signal-to-noise ratio (<x) greater 
than 3. The background is locally estimated from a ring shaped 
region surrounding the star. 

The output catalogs contain for each source an identification 
number, right ascension and declination in degree at J2000, in- 
strumental magnitudes and relative errors, and two SExtractor 
morphological parameters, namely the extraction flag (FLAGS) 
and the isophotal area above the threshold (ISOAREA JMAGE). 
These last parameters were used to clean the catalogs from spu- 
rious detections, such as cosmic ray hits, bad pixels, saturated or 
truncated sources too close to the image boundaries, etc. 

In Figure [3] we show as an example the internal photomet- 
ric errors of all the point-like sources detected in our long-time 
exposures plotted against the magnitude for the RZ filters; the 
relative exponential fits are over-plotted. Table[3]summarises the 
saturation limit and the limiting magnitude achieved at the 10cr, 
5<r and 3<x level in each filter for both the long time exposures 
and the short ones. 

Our final catalog contains celestial coordinates (R.A. and 
Dec. at J2000) and RZ photometry for each point-like object de- 
tected in our survey above the 3<x level. Each entry in the catalog 



3.5. Completeness 

The completeness of our catalogue was estimated in the standard 
way by inserting artificial stars into the images and recovering 
them using the same extraction procedure as for the real objects 
(Sect. [3~4t ; the fraction of recovered artificial objects provides a 
measure of the completeness. 

We used the DAOPHO T II standalone package to perform 
the exercise (IStetsonlll987h . We first use the PSF task to extract 
from the R and Z mosaics the relative PSF models. Then we 
inserted in each mosaic 1000 artificial sources using the addstar 
task; this number should not alter the crowding statistics in the 
images significantly. The profile for the artificial sources was 
generated by using the relative PSF model, while their positions 
are randomly distributed over the entire area of the mosaic and 
their magnitudes range uniformly between the detection and the 
saturation limits in the relative filter (see Table[3]l. 

Figure|4]shows the fraction of recovered artificial objects as a 
function of magnitude from both the deep and shallow mosaics 
for the R and Z filters. The corresponding magnitude limits at 
95% completeness level (C=95%) are reported in Table [3] 

In the absence of extinction (A y-0), assuming a t ypical age 



for Serpens mem bers of 2-6 Myr (Oliveira et al.ll2009l) . a dis- 
tance of 260 pc (IStraizvs et al ]' ll996h and using the theoreti- 
cal isochrones and evolut ionary tracks by iBaraffe et alJ fl998) 
and IChabrier et alJ d2000t) and the completeness limits from 
TableS our optical survey would be complete down to 0.04 M Q 
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Table 3. Number of stellar sources (Nj) detected in both filters above the 3cr level and limiting magnitudes in each filter at 10cr, 
5<r, 3cr and 95% completeness level. 







Filter 


Mag Sat. 


Mag lOcr 


Mag 5cr 


Mag 3cr 


Mag (C=95%) 


Long Time Exposures 


26524 


R 


16.00 


22.90 


23.87 


24.49 


22.30 






Z 


14.50 


21.12 


22.11 


22.73 


19.20 


Short Time Exposures 


3709 


R 


12.50 


19.82 


20.84 


21.50 


18.30 






Z 


10.20 


18.61 


19.41 


19.91 


16.10 



1.2 1 r 



1.0 



0.8 



K 0.6 



0.4 



0.2 



0.0 




20 

Magnitude 

Fig. 4. Completeness (C) plot for extraction of artificial stars 
from our "deep" and "shallow" mosaics for the R and Z bands. 



at the 95% level, i.e. well below the Hydrogen burning limit 
(~0.08 M ). However, as shown bv lHarvev et alJ (l2007l) in their 
Fig. 2, the visual extinction toward the area observed in Serpens 
has a typical value of Ay ~7 mag. Adopting this extinction, we 
estimate a 95% completeness level down to M xQ.l M , 

The full catalog is downloadable from http://cdsarc.u- 
strasbg.fr/cgi-bin/VizieR . 



4. On the disk-less population in Serpens 

The observations presented in this paper are part of the c2d 
complementary work in Serpens. The optical/NIR catalog de- 
scribed in Sect. I3.4l has been merged with the existing X-ray to 
millimeter w avelengths observations collected by the c2d Team 
for Se rpens JGoodmanl 12004 lHarvev et alJl2007t lEnoch et alJ 
l2007t lOliveira et alJ 120091 1201 (C in order to construct com- 
plete spectral energy distributions (SEDs) of the YSOs in the 
observed fields. This dataset allows to anchor disk properties 
of the Serpens YSO sample to the mass, age, and evolution- 
ary status of the central object, which is a critical point in the 
studies of envelopes/d i sks formation and evolution (see, e.g., 
iKundurfhv etail 120061: iMeverl 120091: iPascucci et ai] 120091) . A 
number of papers based on the use of the c2d extensive dataset 
for Serpens and presenting results on disk evolution are in prepa- 
ration. 

As mentioned in Sect.Q] C omeron et alJ(l2009l) have recently 
identified a large population of young members of the Lupus 
dark cloud complex which seems to have lost their inner disks 



on a timescale of a few Myr or less. This discovery poses the 
important question of whether the existence of such disk-less 
young stars is the outcome of specific star-forming conditions in 
Lupus or similar populations exist in other regions. 

XMM-Newton data obtained in April 2007 and April 2008 
of the same field observed by the c2d survey in Serpens revealed 
a sample of new sources, half of which have no c ounterpart in 
the c2 d catalog and are mostly candidate WTTSs dBrown et alJ 
2010). For ~44% of the objects in this sample, we find optical 
and NIR colours consistent with those of young objects with no 
prominent IR excess. Indeed, they occupy in the Z vs. R - Z 
(Fig. [5] upper panel) the same loc us as the YSOs identified by 
the c2d survey (Ha rvey et al.l2 007). which corresponds to an age 
between 1 and 10 Myr acco r ding to pre-m ain sequence (PMS) 
isochrones by iBaraffe et ail d 19981) and IChabrier et all (2000) 



and is consistent wi th the ages of the c2d YSOs found by 
lOliveira et alJ d2009l) . The isochrones have been matched to the 
SDSS photomet ric system fo l lowin g the procedure described in 
Appendix B by ISpezzi et alJ (120071) . On the J - H vs. H - K 
diagram the c2d YSOs follow, as expected, the dwarf stars lo- 
cus with many of them presenting IR excess with respect to this 
locus; the X-ray sources mainly concentrate on the field dwarfs 
locus, as expected for objects with no disk. For the remaining 
objects in the X-ray source sample we could not perform this 
analysis because they are either saturated or not detected in our 
optical/NIR survey, while a few of them are out of the observed 
field. 

To further investigate the disk-less population in Serpens, 
we tried to apply t o our dataset the novel S -parameter method 
by [Comeron et al. (2009), which allows the identification of 
possible young stars and substellar objects based on their op- 
tical/NIR photospheric fluxes, independently of the display of 
signposts of youth, such as IR excess emission or strong Ha 
emission. The S -parameter method simultaneously estimates 
the visual extinction (Ay), effective temperature (T e s) and a 
wavelength-independent scaling factor (S), containing the de- 
pendency on the actual distance and radius of the star, by fit- 
ting a grid of stellar photosphere mode ls to the observed SEP . 
Using the Gala c tic star counts model bv lWainscoat et all dl992l) . 
IComeron et alJ d2009l) demonstrated that members of nearby 
young star forming regions at a given temperature (<4000 K) 
and within a restricted set of ages (<20 Myr) are characterised 
by values of S virtually unreachable by non-members of simi- 
lar temperature, both foreground dwarfs and background gi- 
ants. This method is expected to work reasonably well for the 
detection of cool populations associated to star-forming regions 
located in the ~ 100-300 pc dist ance range from the Sun (see 
Sect. 5 bv IComeron et al . 2009) and, indeed, it has been suc- 
cessfully appli ed to a numb er of nearby star forming regions: 
Lupus (Comer on et al. 2009|). Cha I (Lopez Marti, private com- 
munication), Cr A (Lorjez Marti et al. 2010) and Cha II, where 
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10 My 




4 

R - Z 



2,0 



1.5 



1.0 



0.5 



0.0 



-0.5 




Ay-7 mag 



0.0 



0.5 1.0 
H - K 



1.5 



2.0 



Fig. 5. Z vs. R - Z and J - H vs. H - K diagrams for 
the YSOs in Serpens selected by the c2d survey (circles) 
and the optical/NIR co unter parts of the X-ray sources iden- 
tified by Bro wn et al.l (2010) (squ ares). Upper panel: The 
lines are the PMS isochrones by iBaraffe et alj d!998l) and 
IChabrier etail d2000l) transformed by us into the SDSS photo- 
metric system. The arrow rep resents Ay=2 reddening vector 
dWeingartner & Draind 12001 ). Lower panel: The solid curve 
shows the relation between the colour indices for main sequence 
stars (lower branch) and giants (upper branch), together with the 
relative reddening b ands (dashed fines). The dash-dot line is the 
T Tauri star locus bv lMever et al.l (1 19971) . The object magnitudes 
are derreddened using the typical extincti on in the observed area 
(Ay xT J mag) and the extinction curve bv lWeingartner & D raine 
(2003). 



it provides the same results as in ISpezzi et al.1 (120071) . We ap- 
plied the S -parameter method to the optical/NIR dataset for 
Serpens; we find that in this case the S -parameter histogram 
at any temperature range shows no clear separation between 
the member candidates and the background/foreground contam- 
inants and, as a consequence, the method produces an unreliable 
sample of member candidates, which is also not compatible 
with the XMM-detected sources in the sa me area of the sky 
There are recent indications (see Sect. 2.3 bv lMerin et al . 2008) 
that Galactic models may fail at reproducing stellar counts at 



low galactic latitude and, indeed, Serpens lies very close to the 
galactic plane (b *+5.4 deg). Moreover, Serpens is located at 
260 pc from the Sun, i.e. at the limit of the validity range of 
the S -parameter method. Since all the other clouds on which the 
method has been tested are closer to the Sun and have higher 
galactic latitude than Serpens, we conclude that the unlucky 
combination of distance and position of this cloud prevents the 
S -parameter method from giving useful results. 

5. Conclusions 

We presented an optical/NIR catalog (R and Z filters) of 26524 
point-like sources in 0.96 square degrees in the direction of 
Serpens down to Rw25. These data were collected using the 
WFC camera at the INT within the frame of the Spitzer c2d sur- 
vey in Serpens. The catalog was also complemented with JHKs 
photometry from 2MASS and has been merged with the exist- 
ing X-ray to millimeter wavelengths observations collected by 
the c2d Team for Serpens to study the envelop/disk formation 
and evolution and its dependency on the stellar properties. A 
number of paper based on t his comprehensive ca talog are now 
published or in preparati on dMerin et al ] l2010alrlloiiveiraetalJ 
l2^T^lBTowli^tani20Tol) . 

In this paper we used the optical/NIR catalog to investi- 
gate the disk-less population in Serpens. Because of the dis- 
tance and low galactic latitude o f this cloud, the S -parameters 
method bv lComeron et al.l d2009l) . suitable for the identification 
of young objects independently of any signature of youth, fails 
in distinguishing young cloud members from field contaminants. 
However, a sample of new WTTS candidates has been identified 
in Serpens on the basis of XMM-Newton observations. Our op- 
tical/NIR photometry suggests a very young age (<10 Myr) and 
no NIR excess emission for about 44% of them, supporting their 
WTTS nature. 
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